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! Abstract 

We investigate dynamical chiral symmetry breaking in unquenched QED 3 using the coupled set 
t — \ of Dyson-Schwinger equations for the fermion and photon propagators. For the fermion-photon 

interaction we employ an ansatz which satisfies its Ward-Green-Takahashi identity. We present 
<^ ■ self-consistent analytical solutions in the infrared as well as numerical results for all momenta. In 

Landau gauge, we find a phase transition at a critical number of flavours of Nf rtt ~ 4. In the 
chirally symmetric phase the infrared behaviour of the propagators is described by power laws 
with interrelated exponents. For Nf = 1 and Nf = 2 we find small values for the chiral condensate 
in accordance with bounds from recent lattice calculations. We investigate the Dyson-Schwinger 
equations in other linear covariant gauges as well. A comparison of their solutions to the accordingly 
transformed Landau gauge solutions shows that the quenched solutions are approximately gauge 
covariant, but reveals a significant amount of violation of gauge covariance for the unquenched 



X 

H 

solutions. 
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I. INTRODUCTION 



Over the years, quantum electrodynamics in (2+1) dimensions (QED 3 ) has been studied 
for a variety of reasons. On the one hand it served as a laboratory for investigating nonper- 
turbative phenomena such as dynamical mass generation or confinement in a comparatively 
simple framework devoid of the technical complications of non-abelian gauge theories (for 
reviews see Refs. Qii). O n the other hand QED 3 has regained recent interest due to pos- 
sible applications in condensed matter systems. High-T c cuprate superconductors possess an 
unconventional d-wave symmetry of the pairing condensate. Such a pairing gap has nodes 
at the electronic Fermi surface at which the low energy dispersion becomes linear and thus 
can be described as massless fermions. Since the electronic motion is mainly confined to the 
two-dimensional copper-oxygen planes in these systems an effective low energy description 
of the cuprates in terms of a quantum electrodynamics in two spatial dimensions with two 
massless fermion flavours has been suggested [H, H, E| . In this picture the antiferromagnet- 
ically ordered insulating state of the cuprates would correspond to a state of broken chiral 
symmetry. For this reason there would be considerable interest in a study of the chiral phase 
transition as well as the infrared spectral properties of the fermion propagator in both the 
chirally symmetric and in the ordered phase of QED 3 . 

QED in (2+1) dimensions is a super-renormalisable theory and has an intrinsic mass scale 
given by the dimensionful coupling constant a = Nf e 2 /8. With the help of the photon po- 
larisation n(p) a dimensionless running coupling a = a/ (p[l + U(p)]) = Nf e 2 /(8p[l + II(p)]) 
can be defined which separates the nonperturbative infrared momentum regime from the 
perturbative ultraviolet behaviour 0]. A nonzero fermion mass would provide a second 
mass scale. Various studies of the Dyson-Schwinger equation (DSE) of the fermion propa- 
gator suggest that, in the chiral limit, the interactions generate a dynamical fermion mass 
M(p 2 ) (at least for a small number of fermion flavours), and that this generated mass 
scale M(p 2 = 0) is considerably smaller than the scale defined by the coupling constant a 

It is the smallness of th is g enerated mass scale that poses problems in lattice Monte- 
Carlo simulations of QED 3 0, H H Gij . Finite volume effects are large and the relevant 



signal to determine the chiral phase transition, the dimensionless chiral condensate, is very 
small. Furthermore the presence of an infrared cutoff as such has been shown to reduce the 
value of the critical number of flavours, Nf nt [2oj | . Thus recent studies for the number of 
flavours Nf = 2 [l8| and Nf = 4 jl^, [lj| determined bounds on the chiral condensate, but 
no definite value for Nf Tlt could be extracted. A definite signal for chiral symmetry breaking 
was obtained only for Nf = 1 [19]. Given these problems it seems evident that a continuum 
method is needed to shed light on the infrared properties of QED 3 . 

The DSEs of the propagators of QED 3 have long been investigated employing various 
levels of approximation. Early investigations of the fermion DSE based on a large Nf 
expansion indicated chiral symmetry to be broken only if the number of flavours Nf is 
smaller than a critical value of Nf rit = 32/ V 2 w 3.2 to leading order in Landau gauge [H 
and Nf Tlt = 4/3(32/7r 2 ) including next lo leadingorder corrections in a nonlinear gauge 



2l| . These results have been questioned in Refs. U [1(1 U|, where it was argued that the 



1/iVy-expansion is not an appropriate tool to address these nonperturbative phenomena. 
Using a slightly different truncation of the fermion DSE, it was found that chiral symmetry 
is broken for all values of Nf, although the generated mass scale is exponentially decreasing 



for increasing Nf |ll(. Subsequent work on the coupled DS^s for the fermion and the photon 
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propagator, however, again found chiral symmetry restoration for Nf > Nj , with a value 
of Nf tt between 3 and 4 HI Il5| . All investigations so far are either quenched or employ a 
fermion-photon interaction which manifestly violates gauge symmetry. 

Certainly, gauge invariance is a key property of a local quantum field theory and has to 
play a vital role in these investigations. Reliable results from DSEs can only be expected 
if the fermion-photon vertex respects local gauge symmetry. A necessary (though not suffi- 
cient) condition in this respect is given by the Ward-Green-Takahashi identity (WGTI) (2^ |. 
which determines the 'longitudinal' part of the fermion-photon interaction uniquely in terms 
of the propagator functions (23[ . The remaining transverse part of the vertex is unrestricted 
by the WGTI and has to be determined either directly from the vertex DSE or modelled by a 
suitable ansatz. As the vertex DSE is considerably more complicated to solve than the ones 
for the propagators all work up to now concentrated on the latter strategy. Constraints on 
the structure of the transverse part of the vertex have been derived from gauge covariance 
[24[ HH and multiplicative renormalisability |2(|. Also perturbative vertex corrections put 



constraints on the transverse fermion-photon vertex [27|, |28|, |29|, |3(J, I31 L l32| and p roposals 
for nonperturbative generalisations of these structures have been made |28|, |29|, |30|, |31[ . 

A further important nonperturbative tool to assess the gauge transformation properties 
of a given vertex ansatz are the Landau-Khalatnikov-Fradkin transformations (LKFT) |33| . 
These transformation laws leave the DSEs and the WGTI form invariant and in principle 
allow one to test whether a given ansatz for a vertex is gauge covariant. Such an investiga- 
tion, however, is hampered by the fact that the transformation law for the vertex is quite 
complicated. Therefore an indirect strategy has been applied: one calculates the propagator 
with a given vertex ansatz in the fermion and photon DSE in various gauges and compares 
with the corresponding results from the LKFT of the propag ator 0131113. The 



success of this strategy has been limited by the problem that the LKFT is formulated in 
coordinate space and the necessary Fourier-transform can be carried out analytically only 
for very special cases. 

Our aim in this paper is to make progress in the direction of a gauge covariant solution 
for the propagators of QED 3 . We investigate the coupled set of DSEs for the fermion and 
photon propagator employing a fermion-photon interaction which respects the WGTI. In 
addition, it contains a transverse part that has been shown to respect the LKFT properties 
in the quenched massless case and is thus a good starting point for unquenched QED 3 . With 
the absence of the fermion mass scale in the symmetric phase, the assumption of a power 
law behaviour of the dressing functions in the infrared seems natural. For the same reasons 
similar power laws have been found for the ghost and gluon propagators in quenched and 
unquenched QCD 4 0, H H HJ. We will see how far we can get with this assumption 



here, given the limitations of the chosen truncation scheme. 

The paper is organised as follows: In Sec. |H] we discuss the DSEs for the fermion and 
photon propagator and define our truncation for the fermion-photon interaction. Further- 
more, we recall the ultraviolet behaviour of the propagators as known from methods such 
as the l/iV/-expansion, the loop expansion, and the operator product expansion. In Sec. IHII 
we present an analytical determination of the infrared behaviour of the coupled system of 
fermion and photon equations in the symmetric phase. We show that, within the limits of 
our truncation scheme, the infrared behaviour of the propagators in Landau gauge is given 
by sim p le p ower laws, confirming a longstanding conjecture from perturbative arguments 
|10l IllL |40|. We also calculate the critical number of flavours NJ %t for chiral symmetry 
breaking using these analytic solutions in the infrared, and determine the behaviour of the 
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fermion scalar self-energy close to Nj rit for a simplified version of our fermion-photon vertex. 
Furthermore we investigate the gauge dependence of these power law solutions. 

Numerical results in the broken and symmetric phases are presented in Sec. IIVI In 
passing we re-analyse the quenched fermion and photon DSEs and demonstrate that the 
Curtis-Pennington vertex resolves an inconsistency in determining the chiral condensate 
from the fermion propagator, noted in Ref. We then proceed to solve the unquenched 
system of photon and fermion DSEs employing various ansatze for the fermion-photon ver- 
tex. No further approximations are made. Our results in Landau gauge nicely reproduce 
the analytical results in the ultraviolet as well as in the infrared momentum regime. With 
our most elaborate vertex ansatz the critical number of flavours is NJ' %t w 4. The order 
parameter of the phase transition, the dimensionless chiral condensate (— ( i &ty))/e A is very 
small, of the order of 10~ 3 even in the quenched limit, and decreases exponentially as one 
approaches the phase transition. For N* = 1 and Nf = 2 we find values of the condensate 
in agreement with the lattice bounds 3, We conclude with a discussion of our results 
in Sec. M 



II. THE DYSON SCHWINGER EQUATIONS IN QED 3 



We consider QED 3 with a four-component spinor representation for the Dirac algebra 
and Nf fermions. This allows a definition of chiral symmetry similar to the cases of QED 4 
and QCD 4 . With massless fermions, the Lagrangian has a U(2Nf) "chiral" symmetry, which 
is broken to SU(Nf) x SU(Nf) x £7(1) x U(l) if the fermions become massive 1 . The order 
parameter for this symmetry breaking is the chiral condensate. The question is: is this 
chiral symmetry broken dynamically? We use the set of DSEs to investigate this question. 



A. The fermion and photon propagators 

The DSEs for the photon and fermion propagators in Euclidean space are given by 

d 3 q 



D-J(P) = D'Kp) - Z.Nfe 2 J Tr [ 7|l S(q) T u (q, k) S(k)\ , (1) 

S-\ P ) = S \p) + Z,e 2 / ^3 7/, S(q) T u {q,p) D, v {k) , (2) 

with the momentum routing k^ = q^ — p^. A diagrammatic notation of these equations is 
given in Fig. d 

The general form of the dressed fermion propagator S(p, £) and the photon propagator 
D^(p,0 is given by 

tM(P 2 ,0 + B(p 2 ,Q 



Note that in this formalism a fermion mass term is even under parity. It is also possible to formulate 
QED 3 with two-component spinors; however, in such a formulation there is no "chiral" symmetry, and a 
fermion mass terms breaks parity |4l| . 
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FIG. 1: The Dyson-Schwinger equations of the photon and fermion propagators in diagrammatic 
notation. 

Here £ is the gauge parameter in linear covariant gauges, with £ = denoting Landau gauge. 
The fermion functions A, B, and M depend on the gauge parameter £. On the other hand, 
the vacuum polarisation IT is independent of £. Physical quantities such as the fermion pole 
mass and the chiral condensate are also independent of the gauge parameter. In order to 
keep the notation as clear as possible, we will treat all dependence of the dressing functions 
on £ implicitly from now on. 

The vertex normalisation constant Z\ is related to the fermion wave function normalisa- 
tion Z 2 by a WGTI, Z\ = Z 2 - Since QED 3 is free of ultraviolet divergences, there is no need 
for any renormalisation, though finite renormalisations of the fermion and photon fields are 
possible and leave the physical content of the theory invariant. In our numerical procedure 
we set A(/j, 2 ) = 1 at a large normalisation point fi 2 and determine Z 2 self-consistently. 

B. The fermion-photon vertex 

In general there are several possible strategies to choose an appropriate approximation 
for the fermion-photon vertex T u (q,p) in Eqs. (0) and (J2J). The simplest option would be 
to replace the dressed vertex by the bare vertex j u . However, this violates, among other 
things, gauge invariance and the renormalisation properties of the theory. If one wants to 
preserve these symmetries, one has to use a suitably dressed vertex T u (q,p). 

One way to dress the vertex would be to solve its corresponding DSE. However, the 
vertex DSE contains an unknown four-point function, the fermion-antifermion scattering 
kernel; one has to truncate the infinite set of DSEs somewhere in order to obtain a tractable 
set of equations, and this would only shift the problem up the hierarchy. Furthermore, one 
faces the technical difficulties involved in solving an integral equation in two independent 
momenta, i.e. in three independent variables. 

A different strategy, which we adopt in this paper, is to employ an ansatz for the vertex, 
which has to satisfy at least two requirements: 

(a) it must approach the perturbative form of the vertex for large momenta; 
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(b) it must satisfy the WGTI 



t(q-p)„T u (p,q) = S- 1 (p)-S- 1 (q). (5) 

Condition (a) reflects the fact that QED 3 is an asymptotically free theory as explained in 
the introduction. This condition furthermore implicitly specifies the symmetry properties 
of the vertex, i.e. its behaviour under charge conjugation and Lorenz transformations. 
Condition (b) is dictated by gauge invariance and determines the longitudinal part of the 
vertex. Furthermore, it uniquely fixes the vertex when the two fermion momenta are equal 

. dS-\p) 

r„(p, P ) = ,-^— . (6) 

The two conditions (a) and (b) are necessary but not sufficient to ensure gauge covariance 
of the propagators. We will come back to this point frequently later on. 

Any vertex satisfying condition (b) leads to the following interesting property of the 
fermion equation: With the explicit form, see Eq. P|, of the photon propagator in the 
fermion DSE, the integral on the right hand side can be split into two pieces, 

S-'(P) = S^(p)+Z ie >J ^ 7 „ S{q) pr^Y^pyy - ^) r„( 9 ,„) 

+ Zl ^/(ft 7 " SM ^ r " (? ' p) ' (7) 

with the momentum convention = — p^ for the photon momentum. Due to the 
appearance of the longitudinal projection k u T u (q,p) in the second line of this equation, each 
vertex truncation satisfying the WGTI treats this piece exactly. This will be important later 
on in our infrared analysis. 

A suitable basis to construct a vertex ansatz satisfying the requirements (a) and (b) has 
been given in Ref. [23]. It consists of twelve tensor structures, which can be split up in 
a set of four components, T^ c , completely determined by the WGTI and eight transverse 
components, 

r u (p,q) = r* c (p,q) + r T u (p, q ). (8) 

The WGTI is solved by the Ball-Chiu (BC) construction 



^BC 



Mf) + Atf) . B{P 2 ) ~ B(q 2 ) , 

E^ (p,q) = - + i - 2 - 2 {p + q)i 

2 p z — q z 



The WGTI furthermore constrains the cr^^ + g^) -component of the vertex to be zero. The 
eight transverse components satisfy 

k u T T u (p,q) = 0, T T v (p,p) = Q, (10) 

and are otherwise constrained by condition (a). Much work has been invested to determine 
T^(q,p) in the perturbative region |2?i HE H3, 31, H| to constrain possible 



nonperturbative ansatze, but so far without conclusive results. 
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A minimal ansatz for T^(q,p) ensuring multiplicative renormalisability in four-dimen- 
sional quenched QED has been given by Curtis and Pennington [26| 



v t,cp (t} n) = Mp 2 ) - Mq 2 ) Kp 2 - q 2 h» -{i>-i){p + oU (p 2 + q 2 ) n 1 , 

v {P,q) 2 (p 2 - q 2 ) 2 + (M 2 (p 2 ) + M 2 (q 2 )) 2 ' 1 ' 

Burden and Roberts |24| discovered another favourable property of the Curtis-Pennington 
(CP) vertex T^ p (q,p) = T PC (q,p) + T^' CP (q,p), which holds in quenched massless QED in 
both three and four dimensions: with the help of the LKFT [33j] for the propagators they 
showed that the CP-vertex indeed preserves gauge covariance in these special cases. More 
sophisticated ansatze for the transverse parts of the vertex have been given in Refs. [i^, 



l28l . l29l . l30l . l3l| . Contrary to early expectations, it has been noted [30] that T T (p,q) does 
not vanish in Landau gauge and contains terms that have to be explicitly dependent on the 
gauge parameter £. 



C. The truncated Dyson Schwinger equations 



A numerical investigation of all these ansatze, though highly desirable, is a formidable 
task. Up to now the quenched fermion DSE and the photon DSE of QED^have been studied 
employing a bare fermion photon vertex as well as the BC construction 12 .1241 and the CP- 
vertex j43(. Partly unquenched calculations can be found in Refs. 0,H, id. 111. 44| where 
the photon propagator has been approximated by its 1 /^/-expression. Fully unquenched 
dynamical solutions employing the bare as well as the first part of the BC- vertex have been 
reported in Refs. 0, In our work we will extend these investigations and solve the 
unquenched equations employing the CP construction in the fermion DSE and the BC- 
vertex in the photon equation. The reason for this hybrid choice is the following: The 
transverse term in the CP-vertex has been constructed for quenched QED, i.e. its structure 
is adapted to the kinematical situation in the fermion DSE and it is believed to approximate 
some parts of the real fermion-photon vertex that are important in the fermion DSE. In the 
photon DSE, however, a different kinematical region of the vertex is probed. The T T ' CP -teTm 
leads to divergences here [46| and is thus not a good approximation to the important parts 
of the real vertex in the photon DSE. 

Substituting the CP-vertex into the fermion DSE and taking appropriate traces we arrive 

at 



B{p 2 



Z 2 e 



1 



i + n(p 



(2-Kfk 2 (q 2 A 2 (q 2 ) + B 2 (q 2 )) 
A(p 2 ) + A(q 2 ^ 



2B(q 2 ) + (A(p 2 ) - A(q 2 )) 5(g 2 )Q(p 2 ,g 2 



+ (AA Btf) - ABA(q 2 )) ( + (p 2 + g 2 ) - ( ^ ^ 



2k 2 



A(p 2 ) + A(q 2 



-B(q 2 ) + AA B(q 2 



2\2 



(P ~ g ) 
2k 2 



AB A{q 2 



(P 2 ~ q 2 ) 2 
2k 2 



+ 



q — p 



(12) 
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A{p 2 ) 



Z 2 + Z 2 e 2 



A(q 2 ) 



(2vr) 3 p 2 k 2 (q 2 A 2 (q 2 ) + B 2 (q 2 )) 1 + U(k 2 



A{p 2 ) + A{q 2 ) ( Jf_ {p 2 - q 2 ) 2 
2 + 



2k 2 



+ (A(p 2 ) - A(q 2 )) n(p 2 ,q 2 " 



p z + q l 



AA 



p z + q l 



A- 2 



(P 2 ~ q 2 ) 2 
2k 2 



+ ABM(q 2 ) [-— + (p 2 + q 2 ) 



A{p 2 ) + A(q 2 ) ( p 2 + q 2 (p 2 - q 2 ) 2 



2k 2 



AA- 



(p 2 — q 2 ) ( p 2 — q 2 q 4 — p 



+ 



2k 2 



AB M{q 2 ) 



(p 2 - q 2 ) (q 2 -p 2 ) 2 

2 2k 2 



(13) 



Here we used the abbreviations 

AA 
AB 

n( P 2 ,q 2 ) 



A{p 2 ) - A{q 2 ) 
B(p 2 ) - B(q 2 ) 



p 2 — q 2 



4 4 

p — q 



(p 2 - q 2 ) 2 + (M 2 (p 2 ) + M 2 (q 2 )) 2 



Furthermore we have used the WGTI Z\ = Z 2 for the vertex and wave function renormali- 
sation constants. 

In the photon equation we contract the Lorenz indices with the general tensor 47, 48| 



V<$(p) = 6 



PyPv 
p2 



(14) 



Inserting the BC-vertex in the general expression for the vacuum polarisation, we obtain 

2 2 f d 3 q 1 1 

U{p ) = - 2 e N f J q 2 A 2^ + B 2^ k 2 A 2 {k 2 ) + B 2 {k 2 ) X 

A{q 2 ) + A(k 2 ) ^ fc2)A(g2)A(fc2) + W2{p2j g2? k 2 )B{q 2 )B{k 2^ 

(W 3 (p 2 , q 2 , k 2 )A{q 2 )A{k 2 ) + W 4 (p 2 , q 2 , k 2 )B{q 2 )B{k 2 )) 
(W 5 (p 2 , q 2 , k 2 )A(q 2 )B(k 2 ) + W 6 (p 2 , q 2 , k 2 )B(q 2 )A(k 2 )) \ . 



+ 



+ 



2 




A{q 2 ) - 


A(k 2 ) 


2{q 2 - 


- k 2 ) 


B(q 2 ) - 


B(k 2 ) 


q 2 - 


k 2 



(15) 

The general form of the kernels Wi as well as a discussion of the (weak) dependence of the 
photon equation on the parameter ( of the projector ()14|) can be found in Appendix [XJ 



S 



Here we give the kernels for the special value £ = 3, suggested in Ref. |47( to avoid spurious 
divergences 

3fc 4 , 2 /2 Qq 2 \ 2q 2 3q 4 



Wi(p a ,g 2 l fc a ) = — -k 2 [- 2 + "-\) -1-^ + ^, (16) 



W 2 (p\q\e) = 0, (17) 
^, f, *) = ^ - *« (1 + + *■ (l - %) + - * + * , (18) 

^4 \^ ^4 J y pi J pZ pi 

r* r / 2 2 ,2x - 6fc4 l2 / 4 12 g 2 \ „ V 6g 4 

W 4 (pW,* 2 = ~^ + ^ 2 - + ^- - 2 + -T--T> 19 

p4 y^z ^4 y pZ ^4 

„, , 2 2 , 2 , 3A; 4 )2 / 4 6g 2 \ , 3g 4 

Ws( P 2 , q \k 2 ) = — - k 2 ^ + + 1 + ^ , (20) 
„, , o 2 , 2 , 3A; 4 ,of6q 2 \ , V 3g 4 

H/ 6 p 2 ,g 2 ,fc 2 = — -A; 2 Mf + 1 - JL + J^. 21 
D. l/iV^-expansion and asymptotic behaviour 

Several previous studies of QED 3 have used the 1/iVj-expansion 2 to justify truncations or 
to calculate the asymptotic behaviour of the dressing functions. Our approach does not rely 
on this expansion, but it is interesting to compare our results with the ones obtained in this 
way. We therefore shortly summarise the anticipated behaviour of the dressing functions 
based on the l/iV/-expansion. 

For the vacuum polarisation, one finds for Nf massless fermion flavours to leading order 
in a 1/iVj-expansion [7( 

n p 2 = -t— = - , 22 

8p p 

independently of the value for the gauge parameter £. In the full theory, this expression 
remains valid in the ultraviolet asymptotic limit, as has been demonstrated in quenched 
approximation employing a BC- vertex in the fermion loop ji^. Our numerical results show 
that this is also the case in the unquenched case. In the infrared we will find a modified 
power law for the photon. 

The asymptotic behaviour of the vector self energy to two loop order in quenched approx- 
imation is given by [29] 

if! + e ^ 2 + 3e4 (t_ _ l 

16p 647r 2 p 2 647r 2 p 2 \ 4 3 



^^ 00 ) = 1 + S|_ + _^_ + ^_(^__i) +0 ( a 3) (23) 



Note that to this order the vector dressing function A receives positive corrections in all 
gauges, i.e. A{p 2 ) { 1 for p 2 — ► oo. Furthermore, these corrections to its asymptotic value 
A = 1 behave like inverse powers of the momentum. 



2 This is equivalent to a perturbative expansion for small e 2 while keeping a = Nfe 2 /8 fixed. As QED 3 is 
an asymptotically free theory this expansion will provide correct answers in the ultraviolet. 
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On the other hand, from an unquenched 1 /^/-expansion in Landau gauge, employing a 
bare vertex and the 1/Nf photon propagator given in Eq. (j22j) . one obtains for p < a 



A(p 2 ) = 1 + ln(p/a) , (24) 

up to terms that are regular for p — > 0. Certainly this expression cannot be valid in the 
ultraviolet region (as A(p 2 ) — > 1 for large momenta), nor in the (far) infrared region, re- 
flecting the inconsistency of ignoring vertex corrections. Nevertheless, it has been argued 
Tol 11, 4(J that it could be the first term in the build up of an anomalous dimension 



with 



A (P 2 )=[^)\ (25) 
4 0.135 

^s^ivT^vT' (26) 

in the infrared region. We will come back to this possibility in our infrared analysis in the 
next section. 

Finally, the analysis of the asymptotic behaviour of the scalar dressing function B(p 2 ) 
in the chirally broken phase of massless QED 3 is outlined in Since we are interested 
in dynamical chiral symmetry breaking, which is a purely nonperturbative phenomenon, 
we cannot rely on the 1 /iV/-expansion to obtain an expression for B. Using the operator 
product expansion however, one finds that asymptotically 



*<^oo)=>±lM, (27) 



with A(p 2 — > oo) — ► 1. Thus the chiral condensate (^>^>) can be obtained from the fermion 
propagator in two ways: on the one hand it can be read off from the asymptotic behaviour of 
the B-function and on the other hand it is given by the trace of the propagator in coordinate 
space. In Ref. [HI slight deviations between these two methods have been found. We will 
demonstrate that these deviations do not occur in our truncation. The asymptotic form, 
Eq. (J2ZJ), is reproduced to very good accuracy for a range of values of the gauge parameter 
£ in our numerical analysis. 



III. INFRARED ASYMPTOTIC BEHAVIOUR 



Unlike the situation in a perturbative analysis, where one has a definite starting point 
to work out results order by order, an analysis of the nonperturbative momentum regime 
of QED 3 is based solely on self-consistency and relies therefore on a physically motivated 
ansatz to start with. Inspired by the conjecture of the previous subsection and guided by our 
numerical analysis, our working hypothesis will be that at least in Landau gauge the infrared 
behaviour of QED 3 in the symmetric phase is given by power laws. We will investigate this 
assumption employing different vertex truncations and see how far we can get. Finally we 
will investigate whether our results are gauge covariant. 
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A. Infrared analysis in Landau gauge 



1. The symmetric phase 

For the following analysis we will not use the full CP or BC vertex constructions but only 
the term proportional to 7 M , denoted 1BC, 

This choice has the advantage that the equations are simplified significantly and it already 
contains all qualitative features of the solution employing the full CP/BC- vertex in the 
infrared region, as will be demonstrated by our numerical calculations given in Sec. II VI 
Furthermore, it has the merit that we can solve the DSEs in the infrared analytically. This 
vertex has been considered before in Ref. |15|, though no attempt was made there to solve 
the DSEs analytically. 

The starting point of our investigation is a power law ansatz for the vector dressing 
function 

A{p 2 ) = cp 2 \ (29) 

with the constant c and the power k to be determined self-consistently. We expect this 
ansatz to be valid in the infrared, i.e. in the momentum region p <C a; for p > a the 
function A(p 2 ) rapidly approaches its free form, A(p 2 ) = 1. The integrals on the right hand 
side of the DSEs are dominated by the infrared contributions, coming from the region p < a. 
Thus one can safely substitute the power law for A{p 2 ) and cut off the integrals at p — a. 
Alternatively, one can substitute the power law over the entire momentum range in the 
integrals, provided one keeps track of possible ultraviolet divergences. After integration, the 
resulting power behaviour on the right hand side of the equations then has to match the 
power law on the left hand side. 

Given the ansatz, Eq. (J29|) . we first have to derive the corresponding power law of the 
photon polarisation. After substituting the ansatz (J29j) into the right hand side of the photon 
equation, Eq. (|15j) . the integral can be carried out with the help of Eq. (jB6|) . We arrive at 

n (J) - 7 a 4 r(3/2-«)r(l/2 + «) X _ 2K 

UlBcip } " 22 c n r(3-«)r(i + «) p ' (30) 

= : Z 2 -w 1B cP , 

c 

where we have introduced the dimensionless function Wibc( k ) 

4 r(3/2-«)r(i/2 + «) 
WlBc{K) ='■ n r( 3 -«)r(i + «) • (31) 

For the analysis of the fermion DSE, Eq. ([13)1 . we assume that k > —1/2, as k < —1/2 
only admits the trivial solution A = 1, cf. Sec. 1111 Bl below. With pCa the photon dressing 
is given by 

' ( p l+2 \ (32) 



1 + Tl 1BC (p 2 ) awi B c 
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Together with the power law Eq. ()29|) we then obtain 



c**-Z+ ° fd^a { ^ ^ + q " i- k2 + ( p2 - g2 H \ (33) 

The treatment of this type of equation has been discussed in detail in Refs. 0, HJ for the 
system of ghost and gluon DSEs in QCD 4 . To proceed one has to distinguish two cases, k < 
and k > 0. In the first case the left hand side of the equation becomes singular for p 2 — > 
and has to be matched by a corresponding singularity in the integral on the right hand side. 
In this case the constant term, Z 2 , stemming from the bare propagator is suppressed and can 
simply be discarded. (This case is analogous to the gluon equation in QCD 4 ) . On the other 
hand, if k > the left hand side goes to zero. The renormalisation constant Z 2 thus has to 
be cancelled by a constant term generated by the integral. A straightforward way to deal 
with this situation is to discard the constant term Z 2 and at the same time to eliminate the 
constant term hiding in the integral by employing dimensional regularisation. Furthermore 
we have to eliminate a spurious divergence introduced by employing the power law ansatz 
over the whole momentum range. We are then left with 

2 « = p» ( i + _^ r(,pr(i - «) ) (34) 

wibc Nf vr 2 V2k(1 - 2k) (3 + 2ft) T(3/2 - k)T(1/2 + k) J y ' 

Note that the normalisation factor Z 2 as well as the coefficient c of the power law has been 
dropped out of the equation as expected. The powers of momentum match on both sides 
of the equation thus confirming that the power law is indeed a self-consistent solution of 
the DSEs in the chirally symmetric phase. Equations ()31|) and ()34|) together determine 
the exponent k and therefore completely describe the behaviour of the photon and fermion 
propagators in the infrared in the given truncation scheme. 

From our analysis we find a possible explanation, why the authors of Refs. 

HEQQ did 

not find a phase transition in their truncation scheme: as the feedback from the function A 
onto the vacuum polarisation is not taken into account in their approach, i.e. H{p 2 ) ~ 1/p, 
the right hand side of the DSE for A is proportional to p°, which only matches the left hand 
side iff the A-function becomes a (trivial) constant in the infrared. Thus there is no self- 
consistent power law solution in this truncation scheme. This feedback was first considered 
in Refs. HE Q, H 

An explicit numerical solution of the Eqs. (j31|) and ()34|) is shown in Fig. |21 For the sake 
of comparison we also display the solution for the case of a bare fermion-photon vertex, 
which can be obtained from a similar analysis. Both results are very well fitted by a series 
of powers of 1/Nf. 

0.135 0.090 ... 

Hare = + + 0(l/JVj) , (35) 

0.115 0.044 

KiBC = — + -jys- + 0(l/N}) , (36) 

J f 

which suggests a connection to the l/iVy-expansion. Comparing the first term of our result 
for the bare vertex with Eq. (|25J) we find that the 1/Nf- result is indeed the first term in 
the build up of an anomalous dimension. The additional 1/iVj-term in our fit indicates that 
also loop corrections to the next order in a 1/iVj-expansion sum up and contribute to the 
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0.25 



0.2- 




FIG. 2: Here we display the anomalous dimension of the fermion vector dressing function obtained 
from our infrared analysis, n(Nt), compared with the conjecture from perturbation theory, rj(Nf). 



anomalous dimension. However, one should keep in mind that our calculation is not a 1/Nf 
expansion of the DSEs. It is therefore not surprising that our result for the order 1/Nj 
contribution to the anomalous dimension deviates from that obtained in a 1/Nf expansion 



8(32 - 3tt 2 
0.135 0.022 



V»f = ^W~ -^atT- + 0(1/*?) (37) 



3n 2 N f 9ir 4 N- f 



N f Nj 



0(V*/) • (38) 



Furthermore, the vertex dressing of the lBC-vertex modifies n(Nf) to quite some extent. 

The infrared powers presented in Fig. |21 are not the only solutions of the Eqs. (|3*T|) and 
(134)1 . In the range — 1/2 < k < 1 we found a second solution which is excellently fitted by 

n „ 0.050 0.006 0.028 
"me = O' 5 -^-^-^- (3 9 ) 

However, contrary to the solution ()36j). this solution does depend very heavily on the pro- 
jection method in the photon equation, i.e. it depends on the parameter ( introduced in 
Eq. 1)14)1 . Furthermore, it does not connect to the ultraviolet behaviour of the dressing func- 
tions, i.e. we do not find numerical solutions of the DSEs interpolating between the infrared 
behaviour, Eq. (|39|). and the ultraviolet asymptotic behaviour given in subsection III Dl We 
therefore discard the solution (J39)) in the following. 

We have shown so far that in the symmetric phase the power law ([29)1 leads to a self- 
consistent solution of the fermion and photon DSEs, assuming that the vertex is dominated 
by its 7 M -part. The crucial question is, of course, whether the power law survives when 
additional structure of the vertex is taken into account. That this is indeed the case for 
our choice of the fermion-photon vertex can be shown by a simple dimensional analysis. 
Plugging the power law into the CP- vertex, Eqs. Q and (JUJ), we find that with a vanishing 
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B-function all terms in the vertex depend on combinations of momenta which are of the 
same order p 2n as the leading term. Thus after integration the CP-vertex will only change 
the coefficients of the right hand sides of the fermion and photon DSEs, Eqs. (|3U|) and (|34jh 
but not the general power law behaviour. We thus expect a modified function K CP (Nf) as 
compared to KiBc{Nf) and Kbare{Nf). That these modifications are small will be confirmed 
by our numerical analysis in Sec. II VI below. Further modifications have to be expected 
from including other transverse parts of the fermion-photon vertex and it is by no means 
excluded that k finally becomes negative. We will further discuss this possibility later on in 
Sec. HUB 21 



crit 



2. The chirally broken phase close to 

Next we investigate the chirally broken phase close to the critical value NJ %t of the phase 
transition (assuming for now that there is a chirally broken phase, and a critical value of Nf). 
In this region the dynamically generated fermion mass will be extremely small compared to 
a. The momentum range B(0) <C p <C a will dominate the integral on the right hand side 
of the DSE for the B-function, Eq. ([12)1 . and therefore the chirally symmetric solutions for 
the photon polarisation and the dressing function A 

A(p 2 ) = cp 2K (40) 

rV) = z 2 - c w{k) p - 1 - 2 \ (4i) 

determined in the last subsection, can be substituted in the integral on the right hand side 
for all momenta. Choosing again the case of the lBC-vertex for simplicity we obtain 

B(r? s = 7 f « 1 2B(q) c(p 2 « + q 2 *) 

KP ' 2 N f 7T 3 J q k 2 + Z 2 aWlB c c{K) k 1 ~ 2K c 2 q 2+4K + B(q) 2 2 1 J 

The angular integrals are easily performed, and we obtain 

2 2Z 2 ca 1 7 B(q)(p^ + q 2 ^ ( Z 2 ^f^ + (p + qf^ \ 
[P ) pN f 7r 2 1 + 2k J qq c 2 q 2 +^ + B 2 {q) I z 2 aWlBc(K) + \p - 1 ' 



o 



Q\ 



As we are interested in the momenta p, q <C a the scale a cuts off the integral and the 
logarithm can safely be expanded. Furthermore close to the phase transition the interesting 
region is the one with B 2 (q) c 2 q 2+AK , therefore the equation can be linearised. Thus we 
arrive at 

a 

B(P) = 7~Ylv 2 / dq B(q)iP Z + ^ [m^(p, q))^ . (44) 



The analysis of this type of equation is well known 0, • The integral equation can be 
solved directly by substituting the power law B(p) ~ p b and comparing coefficients on both 
sides. On the other hand, by converting the (nonlinear) integral equation into a differential 
equation one obtains the boundary condition 



P^ + B(p) 



= 0, (45) 

p=a 
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which has to be satisfied by the power law solution. Therefore a nontrivial exponent b has 
to be either complex or b = — 1. 

From the integral equation, Eq. (jSJ), we obtain both the chirally symmetric solution 
B{jp) = and the nontrivial solution 

b = -- + 2k 
2 



± 



'1-4k + 8k 2 ) 



16(1 - 2k) 



-4 




2k) 



16k 2 a 



lu(k) 



2k) 16(1 



2k) 



u(k)Y 



(46) 



where we have used the abbreviation uj(k) := wiBci^^N f. Above a critical value N^ nt 
both solutions of the exponent b are in the interval — 1 < b < and not compatible with the 
boundary condition Eq. ()45|) : The system is in the symmetric phase. Setting the discriminant 
of the outer root in Eq. ()46|) equal to zero and using Eq. ()30|) we find a critical number of 
flavours of 

N J^bc -3.56, (47) 

for the case of the lBC-vertex. For the sake of comparison we also give the result in the 
bare vertex truncation 

3.96. (48) 



jycrit, bare 



Note that for k = we have A = 1, H(p 2 ) = a/» and W\bc = 1; an d consequently recover 
the well-known result from the l/A^j-expansion [8J] 




(49) 



which gives 

Nj UA/Nf =32/vr 2 ^ 3.24. (50) 

(This limit also serves to determine the correct sign in front of the inner root of Eq. (J46j) .) 

Although the critical number of flavours is not too far away from the old 1/Nf result 
there is a clear qualitative difference between solutions from the full (coupled) set of DSEs 
and the one from an l/A^/-expansion: The nonperturbative nature of the DSEs manifests 
itself in the power law solution of the vector dressing function A(p 2 ), i.e. in k ^ 0. Such a 
behaviour can never be obtained in a perturbative expansion and is in marked contrast to 
the assumption A = 1 employed in the leading order l/A^/-expansion. Similar criticism has 



been raised in Refs. 0,0, U, 3|, and is at the origin of a longstanding controversy. 



Having determined the location of the phase transition we now investigate the behaviour 
of the B-function for Nf 



pjcrit f rom b e icw. Abbreviating 
b=-^ + 2K±^f 1BC (N f ,K) 



(51) 



the oscillatory solution of the linearised equation, Eq. 



can be written in the form 



B(p) 



-f 1BC (N f ,K) (ln(p/B(0)) + 6) 



(52) 
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with a phase 5 and the relevant scale for mass generation, B(0), in the logarithm. Plugging 
this solution into the boundary condition Eq. (J43J) we take the limit Nf — > NJ %t and arrive 
at the condition 

l -h BC {N f , k) (ln(p/S(0)) + 5) = nn - -l^-f^Nf, «) . (53) 
It can be shown 0, H| that the value n — 1 gives the lowest vacuum energy. Therefore 

B(0) = ae(w +J ) exp 

We find an exponentially decreasing mass M(0) = B(0)/A(0) at zero momentum close to the 
chiral phase transition. This is in agreement with the numerical findings of Ref . , where 
the same IBC- vertex was used in both the fermion and the photon DSE. A similar expression 
describes the B- function in the bare vertex truncation and in the l/iV/-expansion [71. l8l l21|. 
Our numerical study in Sec. II VI will demonstrate that an exponential decrease of B(0) with 
a modified function f(Nf, k) will also emerge when the CP- vertex is employed. This type of 
exponential behaviour near Nj rit , which is different from the usual first or second order phase 
transition, is reminiscent of a conformal phase transition j52jj]. Strictly speaking however, it 
is not a conformal phase transition, because QED 3 is a super-renormalisable theory with a 
dimensionful coupling constant, and the conformal symmetry is broken in both the chirally 
symmetric and the broken phase. 

B. Infrared analysis in general linear covariant gauges 

Having determined the infrared behaviour of the photon polarisation and the fermion 
dressing functions in Landau gauge we now turn to general linear covariant gauges. In 
the following we will investigate whether our ansatz for the fermion-photon interaction is 
sophisticated enough to generate gauge covariance of the photon polarisation and the fermion 
propagator in the symmetric phase. To this end we will follow the strategy to first re- 
analyse the DSEs in general linear covariant gauges and then compare our findings with 
the corresponding ones from performing a LKFT of our Landau gauge solutions. We will 
investigate whether our truncation allows for power law solutions in general linear covariant 
gauges and whether the LKFT is consistent with such a scenario. 

1 . The photon and fermion DSEs 

In QED the photon polarisation is a gauge invariant object. This is evident from its 
LKFT, given below in Eq. ([63)1 . Thus in general a subtle interplay of the fermion propagators 
and the fermion-photon interaction has to guarantee the invariance of the photon polarisation 
in its DSE. On the perturbative one-loop level, it has been shown recently [3(|, that at least 
in the symmetric phase there have to be terms in the transverse part of the vertex that 
are explicitly dependent on the gauge parameter £. These terms are not constraint by the 
WGTI, and are missing in the vertex truncation investigated in this work. For general 
gauges these terms will be important in the photon DSE and we therefore cannot expect 
the photon polarisation to be gauge invariant at our level of truncation. 



-2tt 



flBc(Nf,K) 



(54) 



16 



Assuming a ^-dependent power law for the vector dressing function of the fermions 



A(p 2 ) = c(0p 



(55) 



and employing the BC-vertex in the photon DSE we end up with the same expression for 
Il(p 2 ) as in Landau gauge, 



n(p 2 ) 



a 



40 



(wibc(k(£)) +W2Bc(k(0)) V 



-l-2«(0 



(56) 



where w\bc has been given in Eq. and w 2 bc abbreviates contributions from the remain- 
ing term of the BC-vertex. The detailed form of W2bc can be calculated analytically, but 
will not be given here as it is not important in the following. We note that on this level of 
truncation gauge invariance of the photon polarisation in the infrared requires the exponent 
k to be independent of the gauge parameter £. It is somewhat surprising that as a result, 
also the infrared behaviour of the vector fermion dressing function A is gauge independent, 
since it is governed by the same exponent k. This may or may not be an inconsistency in 
our truncation as this function A is in general a gauge dependent object. At this stage of the 
investigation one may hope that the dependence of k on £ is weak, leading to qualitatively 
similar results at least in the vicinity of Landau gauge. 

Next we analyse the fermion DSE, Eq. (J7|), with the CP- vertex. In the £-part of the 
equation the vertex is replaced by the inverse fermion propagator according to the WGTI. 
We then obtain 



A(p 2 ) = Z 2 + 



+ 



Z 2 a 
N f n 3 



p2 fc2 ^1 



d 3 q< 



A{p 2 



1 



Tl(k 2 ) p 2 k 2 q 2 A(q 2 ) 



A(p 2 ) + A(q 2 ) + A(p 2 ) 



2 



p 2 q 2 



4 2 
P p 

2P + T 



2k 2 , 

Mi 2 ) p 2 + q 2 ^ 



p A — q z 
p 2 q 2 



2k 2 



2k 2 2 



x 



(57) 



As has been noted in Ref. [2J] the CP-vertex leads to gauge covariant DSEs in the quenched 
massless case. Indeed in this limit, U(k 2 ) = 0, and we have the gauge covariant solution 



1 



A(p 2 



8 7T p 



arctan 



Qlip 



(58) 



of Eq. (|57jl. In Landau gauge, this solution reduces to the trivial solution A{p 2 ) = 1. The 
question is: is there a self-consistent power law solution for the unquenched case? 

Substituting the power laws, Eqs. (|55j) and (|56|). in Eq. (|57|). we have to distinguish 
several cases: 



(a) k > -1/2: 

In this case the photon propagator contributes as 

1 c 



,1+2k 



l + U 1BC (k 2 ) 



Z 2 aw 1B c 



(59) 
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and with substituted power laws the integration of the right hand side of Eq. ([57}) 
leads to 

= z'{ - ^ )+ i f r ( 3/2- K )r(i/2 + K) | 

A^7T 2 ^2 (WlBC + W 2 BC) V "(1 + «)r(l - «) J 

where h(K) abbreviates a combination of T-functions and hypergeometric functions, 
which need not to be specified here. The longitudinal 1/p-term dominates the right 
hand side for all gauges except Landau gauge and again we do not find a self-consistent 
power law to this level of truncation. The only way to obtain such a solution with 
k > —1/2 is the presence of a ^-dependent transverse term in the vertex cancelling 
the ^-dependent longitudinal piece 3 . Then the precise value of k is again determined 
by the coefficients of the remaining terms. 

(b) k < -1/2, k ^ -1: 

In this case the vacuum polarisation vanishes for k — > and the photon propagator is 
proportional to 1/k 2 in the infrared region. Integrating the right hand side of Eq. (|57j). 
the ^-independent part vanishes, and we obtain 

2K 1 Z 2 a r(3/2-«)r(l/2 + «) 

p " P ^ N f 7T r(i + K)r(i- K ) ' 1 ' 

Note that this remaining ^-dependent expression is exact (cf. the comments below 
Eq. (J7J). Matching the coefficients (naively) gives k = —1/2, but due to the divergence 
of the coefficient on the right hand side this is not a solution. There might exist a 
potential solution with k = — 1, which we will consider separately below. Here we 
conclude that there is no nontrivial self-consistent power law solution with — 1 < k < 
— 1/2 within our truncation. Any possible nontrivial solutions in this range have to be 
generated by other transverse parts of the vertex. For example, a term proportional 
to p 2K after integration with appropriate coefficients would be more singular in the 
infrared than the £-piece and lead to a self-consistent nontrivial solution. 

(c) K = -1 

As in case (b) , the vacuum polarisation vanishes for k — > and thus Eq. (J57j) becomes 
effectively quenched. In this case there is at least one solution, namely Eq. ()58j). 
Certainly, this trivial solution will not survive when further parts of the transverse 
vertex will be taken into account, but it serves well in the following to illustrate an 
important point: interestingly, the infrared behaviour of this solution is 

^ 2) = 3 (S) 2 ^ + 0(P4) ' (62) 

for p e 2 £/(87r) and thus k = —1. Nevertheless the corresponding pure power law 
is not a self-consistent solution, as can be seen from Eq. 1)610: the right hand side 
vanishes for k — — 1. The reason for this behaviour is to be found in the appearance of 
the new scale e 2 £/(87r) introduced by the gauge transformation from Landau gauge to 
general linear covariant gauges. For small values of the gauge parameter this new scale 
divides the momentum range < p < a into two regions and in general we cannot 
expect the pure power law to describe the physics in this whole region. 



Based on a less rigorous analysis of the infrared a similar conjecture has already been made in 



0- 
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In all cases considered so far we do not find a self-consistent power law solution in non- 
Landau gauges. Two possible reasons for this behaviour have been identified. Missing 
transverse parts of the fermion-photon vertex could play an important role in these gauges. 
Furthermore the appearance of the new scale e 2 £/(87r) may invalidate a simple power law 
ansatz in the infrared region. In the next subsection we will investigate, whether one can 
derive additional information from the LKFT. 



2. Landau-Khalatnikov-Fradkin transformation of the Landau gauge solution 

Assuming for the moment that the power law for the vector dressing function is qualita- 
tively correct in Landau gauge, we will now use the LKFT to determine the corresponding 
solutions in other gauges. The transformation laws for the photon and fermion propagators 
are most easily specified in coordinate space and we give the transformation rules for the 
propagators in Euclidean space. The photon propagator D /xl/ (x, A) in general gauges can be 
obtained from its transverse Landau gauge form D^x, 0) by the transformation law 

D^(x, A) = D^ix, 0) + d^A(x) , (63) 

with the arbitrary function A. The corresponding transformation law for the fermion prop- 
agator is 

S(x, A) = S(x, 0) e ( A M- A (°)) e2 . (64) 

These transformation laws leave the DSE and the WGTI form invariant. In linear covariant 
gauges and in general dimension d the function A(x) is given by 

/M p -iq-x 
Wf—< < 65 » 

which leads to the familiar form of the photon propagator in linear covariant gauges 



"" v 2 ) p 2 (i + n( P 2 )) +e ' (66) 



in momentum space with the gauge invariant photon polarisation U(p 2 ). Furthermore for 
QED 3 one obtains the transformation law 

S(x,O = S(x,0)e-^ 2 ^\ (67) 

for the fermion propagator in coordinate space. 

In the symmetric phase of Landau gauge QED 3 we found the power law solution 

^'° ) = 7^' (68) 
which leads to the corresponding expression 

d3 P c -i P -x¥ 1 



S(x,0) = J 

(69) 



(27r) 3 p 2 cp 2K ' 
T(l — 2k) sin(K7r) 1 — 2k 7^ 
C 47T 2 k x 3 ~ 2k 
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in coordinate space (all integration formulae used in this subsection are given in Appendix 
IB]) . Applying the transformation ()67|) we transform the propagator to general gauges and 
perform the inverse Fourier-transform 



s( P , 



T(1-2k) sin(KTr) 1 - 2k 



if) 



cAtt 2 
1 2k 



K 



d xe 



ip-x 



li X i c -xge 2 /(87r) 



3-2k 



X" 



(70) 



P z C COS(K7Tj 



2k 



cos[2/t arctan(p87r/ (£e 2 ))] 
[p 2 + (£e7(87r))T 

sin[(2/t — 1) arctan(p87r/ (£e 2 ))] 

_ [p 2 +(^ 2 /(87T)) 2 r 1/2 p(2K-l) 



(71) 



Note that in performing the LKFT with the infrared power law alone we have implicitly 
assumed that contributions from p > a = Nfe 2 /8 have no significant influence on the 
behaviour of the transformed propagator for p<a. Furthermore note that the LKFT has 
introduced a new scale, £e 2 /(87r). In order to be consistent with the previous assumption 
we have to restrict the gauge parameter to small values, i.e. <C £e 2 /(87r) <C a. We then 
obtain two momentum regions of interest where we can expand our solution: 




COS! K7T 



1-Ak 2 



8tt 



2k+2 



for 
for 



(72) 



As expected this expression smoothly connects to the Landau gauge power law when £ — > 0. 
In all other linear covariant gauges we obtain the Landau gauge power law for momenta 
p ^> Below this scale we find surprising agreement with the LKFT of the trivial 

solution, Eq. (JBTJ, i.e. a power k = — 1. 



3. Self- consistent power law solutions in covariant gauges 

If we take the LKFT result based on the power law solution in Landau gauge, Eq. (j72j). 
at face value and combine it with the information we extracted from the infrared analysis of 
the coupled fermion and photon DSEs, Eqs. (which implicitly defines winc and W2bc), 
(I60|) . and two consistent scenarios of massless QED 3 with power law behaviour in the 
infrared are possible: 

I) Landau gauge behaves differently in the extreme infrared than any other linear co- 
variant gauge. In Landau gauge we have a power law with small positive or negative 
values of the exponent k (dependent on the details of the vertex truncation). In other 
gauges we essentially obtain the free solution with k = — 1 for p *|- and the same 

power law as in the Landau gauge for momenta |^ <Cp<Ca. As the region p -C 
can always be gauged away, all interesting infrared physics is already contained in the 
Landau gauge power law. In the DSE this scenario could be realised by a transverse 
term in the vertex which is proportional to the gauge parameter £ and leads to a term 
of the order 1/p 2 in the infrared after integration of the fermion DSE. A subtle can- 
cellation of terms in the photon equation has to guarantee a gauge invariant photon 
polarisation. 
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II) Landau gauge behaves similar to other linear covariant gauges. This entails that the 
solutions with small positive values for k found in Sec. IIII Al are artifacts of the trun- 
cation scheme, and the true solution is a power law in the infrared with an exponent 
k — — 1 for all values of the gauge parameter £ including Landau gauge. Then the 
gauge dependent scale ^ does not distinguish two momentum regions with different 
behaviour of the propagator. In the DSE this scenario requires a transverse term in 
the vertex which does not explicitly contain the gauge parameter £ and leads to a term 
of the order 1/p 2 in the infrared after integration. Then the dressing function A and 
subsequently the photon polarisation would be gauge invariant in the infrared. The 
fermion and photon DSEs would effectively decouple for momenta p<a. 

Both possibilities are consistent with our findings from the LKFT and only a detailed analysis 
of the transverse parts of the fermion-photon vertex can decide which one is realised in QED 3 . 

IV. NUMERICAL RESULTS 

In the previous section we performed in some detail an analytical determination of the 
infrared behaviour of QED 3 close to and above the phase transition assuming a power law 
behaviour of the fermion vector dressing function. Here we present our numerical solutions 
of the unquenched system of DSEs in both the massive and the massless phase, deferring 
the presentation of our results in quenched QED 3 to appendix O 

A. Unquenched results and phase transition in Landau gauge 

Our results for the fermion mass function M(p 2 ), the wave function renormalisation 
Zj(p 2 ) and the photon polarisation H(p 2 ) in unquenched Landau gauge are shown in Fig. 01 
On the left panel we display results obtained with the first term of the BC-vertex (1BC) 
only. On the right panel we give the results obtained with the CP- vertex in the fermion DSE 
and the BC construction in the photon DSE. All results in this section are obtained with the 
Brown-Pennington projection £ = 3 in the photon equation, other choices lead only to minor 
modifications. As expected, in the ultraviolet all curves follow their respective asymptotic 
limits, given in Eqs. (J27)l and (|22j) . In the infrared we find finite, nonzero dressing functions 
throughout the dynamically broken phase. Furthermore, we can clearly see two distinct 
mass scales: The fermion dressing functions have a kink near p = e 2 and a second kink at 
p ~ M(0). Close to the phase transition, these scales are several orders of magnitude apart, 
which makes lattice simulations extremely difficult. 

Above Ny it the functions turn into power laws thus justifying the basic assumption in 
our IR-analysis of Sec. IIII1 The values of the exponents k determined in the analytical 
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TABLE I: The chiral condensate (calculated via Eq. (|73jl) obtained in the lBC-vertex model and 
employing the CP-vertex in the fermion-DSE and the BC-vertex in the photon-DSE. 
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FIG. 3: Shown are the variation of the mass function, the wave function renormalisation and 
the polarisation with the number of flavours in Landau gauge, £ = 0. On the left hand side we 
employed the first term of the BC-vertex (1BC) in both the photon and the fermion DSEs. On the 
right hand side we display the same functions calculated with the CP-vertex in the fermion DSE 
and the BC-vertex in the photon DSE. The scale is set by choosing e 2 = 1. 
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FIG. 4: The phase transition: Analytical vs. numerical results for the scalar fermion dressing 
function, B(p 2 = 0), as function of Nf. Shown are results for four different truncation schemes: 
The 1 /Ay-expansion of Ref. 8], and three different vertex truncations of the coupled photon and 
fermion DSEs. The scale is set by choosing e 2 = 1. 



calculation are reproduced on the 10% level by the numerical results. In the CP/BC- vertex 
case we find k = 0.0315 at Nf = 4.5, which is much closer to the value obtained with a bare 
vertex, K bare = 0.0343, than to the corresponding value k\bc — 0.0278 obtained with the 
lBC-vertex piece. Furthermore, notice that for these unquenched calculations, Zf(p 2 ) > 1, 
at least for Nf > 2, in contrast to our findings in the quenched case (cf appendix IUj) . 

The dynamical mass generation close to the phase transition is studied in Fig. Shown 
are results in four different truncation schemes: The leading order l/A^-expansion of Ref. [1[ 
employs the perturbative expression, Eq. (J2"2~)) . for the photon polarisation and chooses 
A(p 2 ) = 1 together with a bare fermion photon vertex. This is compared to our results from 
the fully unquenched system of DSEs with three different truncations for the vertex: bare, 
1BC and the CP/BC combination. As can be seen from the figure, the numerical results 
follow nicely the corresponding analytical results from Sec. 1111 Al and are in qualitative 
agreement with the findings of Ref. [15]. Our most sophisticated vertex truncation, the 
CP/BC-combination, tends toward a critical value of the number of flavours of Nf Vlt « 4. 

Finally we list our results for the chiral condensate in two different truncation schemes 
for a range of values of Nf in Table U For the CP- vertex, there is no discrepancy between 
the condensate as obtained from the trace of the fermion propagator 

( **> = -nS(0)] = -4 / fy qiAi[ ^ B2(f) , (73) 

and that extracted from the asymptotic behaviour, Eq. (}2Tj) : for the 1BC- vertex there is 
about 5% to 10% difference between the two methods. Listed are the condensates calculated 
from the trace of the fermion propagator. In accordance with the simple estimate given in 



Ref. 54], we find small condensates well below the phase transition. For Nt = 1 and Nt = 2 



23 




10" 2 ^ 4 

10 10 



FIG. 5: Shown are our results for the dressing functions for three different values of the gauge 
parameter £ and Nf = 1. The scale is set by choosing e 2 = 1. 



recent lattice simulations provide upper bounds of the O(10 _3 )e 4 and O(10 _4 )e 4 , respectively 
lii igj | . These bounds are certainly consistent with our values. Thus the combined evidence 
of the DSEs and the lattice Monte-Carlo simulations indicate the presence of dynamical 
chiral symmetry breaking at Nf = 1 and Nf = 2. 



B. Unquenched results in general linear covariant gauges 



Here we present numerical solutions for unquenched QED 3 in linear covariant gauges in 
the CP/BC- vertex truncation scheme. According to our previous discussion in Sec. IIII Bl and 
the numerical results in the quenched approximation in AppendixO we expect artifacts from 
violating gauge invariance. The results presented in Fig. |3] for the case of Nf = 1 flavours 
indicate that this is indeed the case. The photon polarisation functions clearly depend on 
the gauge parameter £. On a quantitative level this can also be seen from the values of the 
chiral condensate displayed in Table [TTJ Induced by the feedback of the photon propagator 
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FIG. 6: Shown are the variation of the mass function, the wave function renormalisation and the 
polarisation with the number of flavours in Feynman gauge, £ = 1. The scale is set by choosing 
e 2 = l. 



on the fermion, the variation of the condensate with the gauge parameter is much larger for 
the unquenched than for the quenched theory. Only the Landau gauge results are consistent 
with the bounds from lattice simulations. 

As expected from the infrared analysis of Sec. 1111 Bl there is no chirally symmetric self- 
consistent power law solution for £ ^ 0. In fact we did not find any self-consistent solution 
in the symmetric phase for £ ^ 0, and consequently the system stays in the chirally broken 
phase at least for values of Nf as large as 7, which is the largest Nf we have investigated. For 
larger values of Nf the numerical analysis becomes increasingly tedious. The corresponding 
dressing functions in Feynman gauge are presented in Fig. |U1 An interesting difference with 
the Landau gauge solutions is that now < Zf(p 2 ) < 1 on the entire momentum range, in 
contrast to the unquenched Landau gauge solutions, for which Zf(p 2 ) > 1 on part (Nf = 1) 
or all (Nf > 2) of the momentum range. Also note that now there appears to be only one 
scale at which the generated fermion mass function function M(p 2 ) has a kink. 
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gauge parameter 






-(**)/(10 _ 


e ) 








A/* — 


Nt — 1 


N f = 2 


N f = 3 


Nt — 4 


Nt — 5 


Nt — 6 


e = o 


333 


121 


13 


0.026 


? 








f = 0.5 


340 


165 


79 


39 


23 


15 


11 


e = i 


351 


202 


108 


74 


55 


37 


29 


e = 2 


356 


259 


189 


143 


107 


92 


77 



TABLE II: The chiral condensate (calculated via Eq. Q73))) obtained in the CP/BC-vertex trun- 
cation in different gauges. 



V. SUMMARY AND CONCLUSIONS 

In this work we have investigated the chiral phase transition of QED 3 in the Green's 
functions approach. Employing different ansatze for the fermion-photon vertex we have 
solved the coupled set of Dyson-Schwinger equations for the fermion and photon propaga- 
tors. No other approximations have been made in our numerical calculations. In addition, 
the infrared behaviour of the propagators close to the phase transition and in the sym- 
metric phase has been investigated employing methods that have been successfully used 



previously in four-dimensional QCD |36j, 1221 12a |39|, |5Cj and QED |52j. Special care has 



been taken to preserve gauge invariance as much as possible. The chosen vertex ansatz 
satisfies the Ward-Green-Takahashi identity. Furthermore it has the correct properties un- 
der Landau-Khalatnikov-Fradkin transformations in the special case of massless quenched 
QED. Nevertheless, our results indicate that this is not enough: the resulting vacuum po- 
larisation is not gauge invariant, nor is the chiral condensate. Clearly, further structure in 
the transverse part of the vertex is needed to properly ensure gauge covariance. 

In Landau gauge we find a self-consistent power law solution for the photon polarisation 
and the vector fermion dressing function in the infrared region in the symmetric phase. With 
a bare fermion-photon vertex the anomalous dimension k is directly related to the coeffi- 
cients of a well known result from the 1 / AVex pan sion. We thus confirmed a longstanding 
conjecture from the renormalisation group (loL 11, 13- We would like to emphasise, how- 



ever, that such a power law solution is genuinely nonperturbative in nature and cannot be 
obtained to any finite order in perturbation theory or the 1/A^-expansion. The dependence 
of k on the number of flavours is modified by nonperturbative contributions in the vertex 
dressing. We find small positive values of k for all vertex dressings considered so far, in- 
dicating a vanishing propagator function A(p 2 ) for p 2 — > 0. This is not what one would 
expect on physical grounds 4]. It remains to be seen whether further contributions from 
the transverse part of the fermion-photon vertex are capable to drive k to negative values. 

In the chirally broken phase, the power law behaviour gets modified by the (dynamically 
generated) fermion mass, which effectively acts as an infrared cutoff. We have determined the 
dependence of the chiral condensate and of the scalar fermion dressing function B(p 2 = 0) 
on the number of flavours and found an exponential decrease close to the phase transition. 
If this behaviour turns out to be correct in the full theory, one can hardly hope to be able 
to determine the critical number of flavours from lattice Monte-Carlo simulations. A sign 
pointing in this direction is found for small Nf\ The chiral condensate is very small compared 
to the natural mass scale e 2 . Furthermore it agrees with the values recently determined on 
the lattice [18, 19]. The qualitative behaviour of the B-function and the condensate does 
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not depend on our choice for the vertex ansatz, though there are quantitative differences, in 
particular near the critical number of flavours, Nj rit . The results with our most sophisticated 
vertex suggest a critical number of flavours Nj rit « 4. 

In other linear covariant gauges we find a completely different picture. No indications 
for a phase transition are seen in our numerical analysis. The value of the condensate 
becomes heavily dependent on the gauge parameter, and for £ 7^ our results exceed the 
bounds set by lattice simulations. It appears as if no self-consistent power law solutions 
exist in the symmetric phase. A possible explanation for this fact can be found with the 
help of the Landau-Kalatnikov-Fradkin transformation laws. We find that, given a power 
law solution in Landau gauge with exponent k, the gauge transformed propagator has the 
same anomalous dimension only for momenta |j- <C p C a, whereas below this scale it 
has an anomalous dimension k — — 1. Such a solution cannot be found from the Dyson- 
Schwinger equations with our vertex truncation. We thus conclude that further transverse 
structure in the vertex is mandatory to obtain gauge covariant solutions for the propagators 
of QED 3 . These extra terms could allow for three possible scenarios in the symmetric phase 
of QED 3 : (a) terms explicitly proportional to the gauge parameter £ could lead to solutions 
in accordance with the Landau-Khalatnikov-Fradkin transformation of the Landau gauge 
solutions; or (b) they could invalidate all Landau gauge solutions found so far; or (c) they 
could allow for solutions in the symmetric phase that are not power laws in the infrared 
region. Although we do believe that the scenario (a) is most likely to be realised we cannot 
exclude the other possibilities so far. A thorough investigation of all existing proposals for 
the transverse structure of the fermion-photon vertex in this respect seems desirable and 
will be carried out in future work. 
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APPENDIX A: REGULARISATION OF THE PHOTON DSE 

The aim here is to provide a formulation for the photon equation that can be used not 
only in combination with the Brown-Pennington projector, ( = 3 (cf. Eq. (|14p). but with 
general values of (. The dependence of the Landau gauge vacuum polarisation on ( then 
provides a measure of the violation of transversality in the photon equation. 

To see the problems arising with ( ^ 3 we analyse the ultraviolet behaviour of the photon 
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DSE, Eq. (|15j). The kernels Wi are then given by 
W 1 (p 2 ,q 2 ,k 2 ) 
W 2 (p 2 ,q 2 ,k 2 ) 



(k* , .j/i-c 20A , , (i-C)? 2 , Cq 
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(A2) 



(A3) 
(A4) 
(A5) 



W^ 2 ,g 2 ,A; 2 ) = 
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< - l)? 2 , (q 
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Let us first concentrate on the term proportional to the kernel W\. In the ultraviolet we 
can use the angular approximation A(q 2 ) = A(k 2 ) = 1, furthermore the momenta q 2 and k 2 
are larger than all fermion masses. One angular integral yields a trivial factor 2tt and the 
integral is dominated by the part q 2 > p 2 



U uvl (p 2 " 



OO 7T 

J dqq 2 J d6sm(6)-^x 



9 2 N f 

(27T) 2 

V 



P 



P 



" 1 + 



P 



P 



(A7) 



Now we perform the angular integrals according to Eqs. (|B1|) through ()B5|) and expand the 
resulting logarithm for momenta q 2 ^> p 2 . To leading order we obtain 



^iV/ /\ (2(3-C) , (-10-20 
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(/(/ 



3p 2 
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+ 0(p 2 /^ 
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which is convergent iff £ = 3 but contains a linear divergence for all other values. Treating 
all other terms in the fermion loop in the same way we arrive at the expression 
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which contains linear divergences proportional to (3 — Q at various places. In order to 
eliminate these terms we subtract appropriate expressions from the kernels W% given in 



Eqs. through (|2ip. This results in the modified kernels 

W x {p\q\k") = W.ip^q^k*)- 2 ^; , (A10) 

W 2 (p 2 ,q 2 ,k 2 ) = 0, (All) 

W 3 (p 2 , q 2 , k 2 ) = W,(p 2 } q 2 , k 2 ) + 8g2 ^ (3 2 ~ C) , (A12) 

op 

W 4 (p 2 , q 2 , k 2 ) = W 4 (p 2 , q 2 , k 2 ) - %k2{ l~ , (A13) 

W 5 (p 2 , q 2 , k 2 ) = W 5 (p 2 , q 2 , k 2 ) + 4fc2( 3 3 p2 ~ , (A14) 

W,(p 2 , q 2 , k 2 ) = W 6 (p 2 , q 2 , k 2 ) + 4fc2( 3 " . (A15) 



Based on our analytical infrared calculus as well as on our numerical calculations we inves- 
tigated the dependence of the solutions on the projection parameter ( and found very small 
effects not affecting any of our conclusions in the main body of this work. 



APPENDIX B: ANGULAR AND RADIAL INTEGRALS 

In d = 3 dimensions the following angular integrals are needed for the UV-analysis of the 
photon equation 

de n^- = (B1) 

de^ = -\n(p^-), (B2) 



o 







k 2 pq \\p — q\ 
d0sin(0) = 2, (B3) 



o 



7T 



J d6sm(9)k 2 = 2(p 2 + q 2 ), (B4) 
o 

7T 

/2f) 
d9 sin(0) k 4 = 2p 4 + 2g 4 + —p 2 q 2 , (B5) 

o 

where k 2 = (q — p) 2 = p 2 + q 2 — 2pqcos(9). For the IR-analysis of the coupled system we 
need the following integrals 

f M n 1 _ *- d /2 ( ,, d/ 2-a- b r(d/2 - g)T(d/2 - b)T(a + b - d/2) 
J q {q 2 Y{k 2 f ~ [P ' T(a)T(b)T(d - a - b) ' 1 ' 
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and for the Fourier-transforms necessary for the LKFT we need 



d6 sin 6 cos 6 e ±ipxcos 9 = 



cos(px) sin(pa;) 

(px) 2 



px 



dxx b sin (ax) 



oo 



r(i + &) . 



r(i + 6) 



,1+6 



COS 



7T 



7T 



(B7) 



'1 + ^2 > 0<|6 + 1|<1, [3.761(4)] (B8) 



(1 + 6)- , <(& + !)<!, [3.761(9)](B9) 



dxx sin(ax) e 



r(i + 6) 



dxx b cos(ax) e 



(a 2 + c 2 )( 1+b )/ 2 

r(i + 6) 



sin 



■ cos 



'1 + 6) arctan ^— ^ 



'1 + 6) arctan I - 



, 6 > -2, c> , 
[3.944(5)] (BIO) 



, 6 > -1, c> 0, , 



(a 2 + c 2 )( 1+6 )/ 2 
o 

[3.944(6)] (Bll) 

where the numbers in square brackets refer to the corresponding equations in Ref. [Hif . 



APPENDIX C: NUMERICAL RESULTS IN QUENCHED APPROXIMATION 

QED 3 in quenched approximation employing the BC-vertex in the fermion and photon 
equation has been investigated in detail in Refs. [l^, IH1 (for recent work see Ref. |35| 
and references therein) . In the chirally broken phase, the gauge dependence of the photon 
polarisation as well as the chiral condensate was found to be rather weak for the condensate 
and uncomfortably large for the photon polarisation. For the condensate a discrepancy 
between the value extracted from the asymptotic form of the scalar fermion self-energy, see 
Eq. (J27j) . and the value obtained from the trace of the fermion propagator, see Eq. (|73j). has 
been found. All we have to add to these investigations is an answer to this last problem. 
As can be seen from Table 11111 adding the CP term to the BC-vertex in the fermion DSE 
removes this discrepancy and slightly reduces the gauge dependence of the condensate. 



(-(^})/(10- 3 e 4 ) 
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0.5 




= 1 




asympt. 


-Tr[S(0)] 


asympt. 


-Tr[S(0)] 


asympt. 


-Tr[S(0)] 


BC-vertex 


3.34 


3.24 


3.54 


3.25 


3.67 


3.26 


CP-vertex 


3.29 


3.29 


3.34 


3.35 


3.46 


3.46 



TABLE III: The chiral condensate extracted from the asymptotics of the fermion mass function, 
see Eq. (|27|). and obtained by taking the trace of the propagator, for different values of the gauge 
parameter £, all in the quenched (Nf = 0) approximation. The units are 10 -3 e 4 . 
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FIG. 7: Here we display the fermion mass function, M(p 2 ), and the wave function renormalisation, 
Zf(p 2 ) = l/^4(p 2 ) in quenched approximation for three different values of the gauge parameter £. 
The scale is set by choosing e 2 = 1. 




FIG. 8: Here we display the photon polarisation H(p 2 ) as calculated from the photon DSE using 
Nf = 1 and the quenched fermion propagator functions without back-coupling for three different 
values of the gauge parameter £. The scale is set by choosing e 2 = 1. 



The corresponding numerical solutions are displayed in Figs. [7| and |H1 For the fermion 
propagator, no qualitative difference between Landau gauge and gauges with non-vanishing 
gauge parameter is found. Notice that < Zf(p 2 ) < 1 on the entire momentum range, as 
one would expect from quenched perturbative theory, eq.(|23|l. For the photon polarisation 
we find sizable violations of gauge invariance employing the BC-vertex in accordance with 
Ref. (4^ . Adding the CP term also in the photon DSE does not help in this respect and 
furthermore introduces spurious divergences. 
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